GENERATING A LIBRARY OF SIMULATED-DIFFRACTION SIGNALS 
AND HYPOTHETICAL PROFILES OF PERIODIC GRATINGS 

BACKGROUND 

1 . Field of the Invention 

[0001 ] The present invention relates to wafer metrology, and more particularly to 

generating a library of simulated-diffraction signals for use in optical metrology to 
determine the profile of periodic gratings. 

2. Related Art 

[0002] In semiconductor manufacturing, periodic gratings are typically utilized 

for quality assurance. For example, one typical use of such periodic gratings includes 
fabricating a periodic grating in proximity to the operating structure of a semiconductor 
chip. By determining the profile of the periodic grating, the quality of the fabrication 
process utilized to form the periodic grating, and by extension the semiconductor chip 
proximate the periodic grating, can be evaluated. 

[0003] Optical metrology can be utilized to determine the profile of a periodic 

grating. In general optical metrology involves directing an incident beam at the periodic 
grating, and measuring the resulting diffraction beam. The characteristics of the 
measured-diffraction beam (i.e., a measured-diffraction signal) is typically compared to a 
library of pre-determined diffraction signals (i.e., hypothetical-diffraction signals) that are 
associated with known profiles. When a match is made between the measured-diffraction 
signal and one of the hypothetical-diffraction signals, then the profile associated with the 
matching hypothetical-diffraction signal is presumed to represent the profile of the 
periodic grating. 

[0004] In general, the process of generating a hypothetical-diffraction signal 

involves performing a large number of complex calculations, which can be time and 
computationally intensive. The amount of time and computational capability and 
capacity needed to generate hypothetical-diffraction signals can limit the size and 
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resolution (i.e., the number of entries and the increment between entries) of the library of 
hypothetical-diffraction signals that can be generated. 

SUMMARY 

[0005] In one exemplary embodiment, a library of simulated-diffraction signals 

and hypothetical profiles of a periodic grating is generated by generating diffraction 
calculations for a plurality of blocks of hypothetical layers. A diffraction calculation for 
a block of hypothetical layers characterizes, in part, the behavior of a diffraction beam in 
the block of hypothetical layers. Each block of hypothetical layers includes two or more 
hypothetical layers, and each hypothetical layer characterizes a layer within a 
hypothetical profile. The diffraction calculations for the blocks of hypothetical layers are 
stored prior to generating the library. The simulated-diffraction signals to be stored in the 
library are then generated based on the stored diffraction calculations for the blocks of 
hypothetical layers. 

DESCRIPTION OF DRAWING FIGURES 
[0006] The present invention can be best understood by reference to the following 

description taken in conjunction with the accompanying drawing figures, in which like 
parts may be referred to by like numerals: 

Fig. 1 depicts an exemplary optical-metrology system. 
Fig. 2 depicts another exemplary optical-metrology system. 
Figs. 3 A - 3E depict various hypothetical profiles of periodic 

gratings. 

Fig. 4 depicts an exemplary cache and library. 

Figs. 5 A and 5B depict an exemplary process of generating 
simulated-diffraction signals. 

Figs. 6A and 6B depict another exemplary process of generating 
simulated-diffraction signals. 

Figs. 7 A and 7B depict still another exemplary process of 
generating simulated-diffraction signals. 



Fig. 8 depicts an exemplary computer system. 

Fig. 9 shows a section of a periodic grating. 

Figs. 1 OA - IOC depict cross sectional views of a portion of a 
periodic grating. 

Fig. 1 1 depicts an exemplary coupled-wave analysis. 

Fig. 12 depicts another exemplary coupled-wave analysis. 

Fig. 13 depicts an exemplary process of generating simulated- 
diffraction signals. 

Fig. 14 depicts another exemplary process of generating simulated- 
diffraction signals. 

Fig. 15 is another exemplary computer system. 

Figs. 16A - 16D depict exemplary ridge profiles. 

Fig. 17 depicts an exemplary aggregate-layer trapezoid. 

Fig. 18 depicts an exemplary process of generating a library of 
simulated-diffraction signals. 

DETAILED DESCRIPTION 
[0007] The following description sets forth numerous specific configurations, 

parameters, and the like. It should be recognized, however, that such description is not 
intended as a limitation on the scope of the present invention, but is instead provided as a 
description of exemplary embodiments. 

[0008] With reference to Fig. 1, a periodic grating 145 is depicted on a 

semiconductor wafer 140. As depicted in Fig. 1, wafer 140 is disposed on a process plate 
180, which can include a chill plate, a hot plate, a developer module, and the like. 
Alternatively, wafer 140 can also be disposed on a wafer track, in the end chamber of an 
etcher, in an end-station or metrology station, in a chemical mechanical polishing tool, 
and the like. 



[0009] As described earlier, periodic grating 145 can be formed in test areas on 

wafer 140 that are proximate to or within an operating structure formed on wafer 140. 
For example, periodic grating 145 can be formed adjacent to a device formed on wafer 
140. Alternatively, periodic grating 145 can be formed in an area of the device that does 
not interfere with the operation of the device. As will be described in greater detail 
below, the profile of periodic grating 145 is obtained to determine whether periodic 
grating 145, and by extension the operating structure adjacent periodic grating 145, has 
been fabricated according to specifications. 

1 . Optical Metrology 

[001 0] In one exemplary embodiment, optical-metrology can be utilized to obtain 
the profile of periodic grating 145. As depicted in Fig. 1, an optical-metrology system 
100 can include an electromagnetic source 120, such as an ellipsometer, reflectometer, 
and the like. Periodic grating 145 is illuminated by an incident beam 110 from 
electromagnetic source 120. Incident beam 110 is directed onto periodic grating 145 at 
an angle of incidence 9i with respect to normal n of periodic grating 145. Diffraction 
beam 115 leaves at an angle of 9 d with respect to normal n . In one exemplary 
embodiment, the angle of incidence Q } is near the Brewster's angle. However, the angle 
of incidence 0i can vary depending on the application. For example, in an alternative 
embodiment, the angle of incidence 0i is between about 0 and about 40 degrees. In 
another embodiment, the angle of incidence is between about 30 and about 90 degrees. 
In still another embodiment, the angle of incidence is between about 40 and about 75 
degrees. In yet another embodiment, the angle of incidence 9* is between about 50 and 
about 70 degrees. 

[001 1] Electromagnetic source 120 can include focusing optics to control the spot 

size of incident beam 1 10. In one embodiment, spot size of incident beam 1 10 can be 
reduced to less than the size of the test area on wafer 140 that contains periodic grating 
145. For example, a spot size of about 50 \xm by 50 \xm, or smaller, can be used. 
Additionally, electromagnetic source 120 can include a pattern recognition module to 
center the spot in the test area on wafer 140. 



[0012] As depicted in Fig. 1 , diffraction beam 1 1 5 is received by detector 1 70. 

When electromagnetic source 120 is an ellipsometer, the magnitude Q¥) and the phase 
(A) of diffraction beam 115 are received and detected. When electromagnetic source 120 
is a reflectometer, the relative intensity of diffraction beam 1 15 is received and detected. 

[001 3] To determine the profile of periodic grating 145, optical-metrology system 
100 includes a processing module 190, which converts diffraction beam 115 received by 
detector 170 into a diffraction signal (i.e., a measured-diffraction signal). Processing 
module 190 then compares the measured-diffraction signal to the simulated-diffraction 
signals stored in a library 1 85. As will be described in greater detail below, each 
simulated-diffraction signal in library 185 is associated with a theoretical profile. When a 
match is made between the measured-diffraction signal and one of the simulated- 
diffraction signals in library 185, the theoretical profile associated with the matching 
simulated-diffraction signal is presumed to represent the actual profile of periodic grating 
145. The matching simulated-diffraction signal and/or theoretical profile can then be 
utilized to determine whether the periodic grating has been fabricated according to 
specifications. 

[001 4] Although optical-metrology system 100 has been described and depicted 

thus far as having a single electromagnetic source 120 and a single detector 170, it should 
be recognized that optical-metrology system 100 can include any number of 
electromagnetic sources 120 and detectors 170. For example, with reference to Fig. 2, 
optical-metrology system 100 is depicted with two electromagnetic sources 120A and 
120B, and two detectors 170A and 170B. Two incident beams 110A and HOB are 
generated by electromagnetic sources 120 A and 120B, with each incident beam 1 10A 
and HOB including two polarizations of the electromagnetic wave to allow for both 
intensity and phase measurements. 

[001 5] As depicted in Fig. 2, an optical fiber 125 can transport broadband 

radiation from a source (not shown) to the electromagnetic sources 120A and 120B. A 
switching mechanism 200 can direct the broadband radiation alternately to the two optic 
fiber branches 125 A and 125B, which lead to the electromagnetic sources 120 A and 



120B, respectively. In a preferred embodiment, the switching mechanism 200 is a wheel 
with an opening in one semicircle, so that rotation of the switching wheel 200 by 180 
degrees produces a switching of the optic fiber branch 125 A or 125B to which the 
radiation is directed. In this preferred embodiment, the incident beams are directed to 
wafer 140 so that the angles of incidence 9n and 9 i2 are approximately 50 and 70 degrees. 
Other angles of incidence are also contemplated as useful, and angles between zero and 
90 degrees are possible. In the present embodiment, incidence Q\\ and 0i 2 are chosen to be 
disparate enough that redundancy of information is not an issue, while both Q\\ and 0*2 are 
chosen to be close enough to Brewster f s angle that the incident beam is sensitive to the 
grating features. 

[001 6] As depicted in Fig. 2, diffraction beams 1 1 5 A and 1 1 5B leave the wafer at 

angles 9<n and 9d2- Diffraction beams 1 15 A and 1 1 5B are then received by detectors 
170A and 170B. When electromagnetic sources 120A and 120B are ellipsometers, the 
magnitude Q¥) and the phase (A) of diffraction beams 1 15A and 1 15B are received and 
detected. When electromagnetic source 120 is a reflectometer, the relative intensity of 
diffraction beams 1 15A and 1 15B is received and detected. 

[001 7] It should be understood that multiple detectors 170A and 170B are not 

necessary to perform multiple-angle reflectometry or ellipsometry measurements. 
Instead, a single electromagnetic source 120 and detector 170 can perform a first 
measurement for an angle of incidence of On, and then be moved to perform a second 
measurement for an angle of incidence of 0i2. 

2. Library of Simulated-Diffraction Signals and Theoretical Profiles 

[001 8] As described above, with reference again to Fig. 1, library 185 includes 

simulated-diffraction signals that are associated with theoretical profiles of periodic 
grating 145. More particularly, in one exemplary embodiment, library 185 includes 
pairings of simulated-diffraction signals and theoretical profiles of periodic grating 145. 
The simulated-diffraction signal in each pairing includes theoretically generated 
reflectances that characterize the predicted behavior of diffraction beam 115 assuming 
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that the profile of the periodic grating 145 is that of the theoretical profile in the 
simulated-diffraction signal and theoretical profile pairing. 

[001 9] The set of hypothetical profiles stored in library 1 85 can be generated by 

characterizing a hypothetical profile using a set of parameters, then varying the set of 
parameters to generate hypothetical profiles of varying shapes and dimensions. The 
process of characterizing a profile using a set of parameters can also be referred to as 
parameterizing. 

[0020] For example, as depicted in Fig. 3 A, assume that hypothetical profile 300 

can be characterized by parameters hi and wl that define its height and width, 
respectively. As depicted in Figs. 3B to 3E, additional shapes and features of 
hypothetical profile 300 can be characterized by increasing the number of parameters. 
For example, as depicted in Fig. 3B, hypothetical profile 300 can be characterized by 
parameters hi, wl, and w2 that define its height, bottom width, and top width, 
respectively. Note that the width of hypothetical profile 300 can be referred to as the 
critical dimension (CD). For example, in Fig. 3B, parameter wl and w2 can be described 
as defining the bottom CD and top CD, respectively, of hypothetical profile 300. 

[0021] As described above, the set of hypothetical profiles stored in library 185 

(Fig. 1) can be generated by varying the parameters that characterize the hypothetical 
profile. For example, with reference to Fig. 3B, by varying parameters hi, wl, and w2, 
hypothetical profiles of varying shapes and dimensions can be generated. Note that one, 
two, or all three parameters can be varied relative to one another. 

[0022] With reference again to Fig. 1, the number of hypothetical profiles in the 

set of hypothetical profiles stored in library 185 depends, in part, on the range over which 
the set of parameters and the increment at which the set of parameters are varied. The 
range and increment can be selected based on familiarity with the fabrication process for 
periodic grating 145 and what the range of variance is likely to be. The range and 
increment can also be selected based on empirical measures, such as measurements using 
AFM, X-SEM, and the like. For a more detailed description of determining the number 
of parameters to use in characterizing a hypothetical profile and the range over which to 



vary the parameters, see U.S. Patent Application Ser. No. 09/907,488 titled 
GENERATION OF A LIBRARY OF PERIODIC GRATING DIFFRACTION 
SIGNALS, filed on July 16, 2001, the entire content of which is incorporated herein by 
reference. 

[0023] As described above, the simulated-diffraction signals stored in library 185 

can be generated by generating a simulated-diffraction signal for each hypothetical 
profile in the set of hypothetical profiles stored in library 185. As will be illustrated in 
more detail below, the process of generating a simulated-diffraction signal can involve 
performing a large number of complex calculations. Additionally, as the complexity of 
the theoretical profile increases, so does the number and complexity of the calculations 
needed to generate the simulated-diffraction signal for the theoretical profile. 
Furthermore, as the number of theoretical profiles increases, so does the amount of time 
and processing capacity and capability needed to generate library 185. 

[0024] As such, in one exemplary embodiment, a portion of the calculations 

performed in generating simulated-diffraction signals can be stored as intermediate 
calculations prior to generating library 185. The simulated-diffraction signals for each 
hypothetical profile in library 185 can then be generated using these intermediate 
calculations rather than performing all of the calculations needed to generate a simulated- 
diffraction signal for each theoretical profile in library 185. 

[0025] More particularly, in one exemplary embodiment, diffraction calculations 

are generated for a plurality of blocks of hypothetical layers. A diffraction calculation for 
a block of hypothetical layers characterizes, in part, the behavior of diffraction beam 115 
in the block of hypothetical layers. Each block of hypothetical layers includes two or 
more hypothetical layers, and each hypothetical layer characterizes a layer within a 
hypothetical profile. The diffraction calculations for the blocks of hypothetical layers are 
stored prior to generating library 185. The simulated-diffraction signals to be stored in 
library 185 are generated based on the stored diffraction calculations for the blocks of 
hypothetical layers. 
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[0026] For example, with reference to Fig. 4, a plurality of hypothetical layers 

402 (e.g., hypothetical layer 402.1 - 402.12) can be grouped together as a block of 
hypothetical layers 404 (e.g., block of hypothetical layers 404.1 - 404.4). Each 
hypothetical layer 402 characterizes a layer of a hypothetical profile. 

[0027] In the present example, hypothetical layers 402.1, 402.2, and 402.3 are 

depicted as being grouped together to form block of hypothetical layers 404.1. 
Hypothetical layers 402.4, 402.5, and 402.6 are depicted as being grouped together to 
form block of hypothetical layers 404.2. Hypothetical layers 402.7, 402.8, and 402.9 are 
depicted as being grouped together to form block of hypothetical layers 404.3. 
Hypothetical layers 402.10, 402.1 1, and 402.12 are depicted as being grouped together to 
form block of hypothetical layers 404.4. In this manner, a number of blocks of 
hypothetical layers 404 of various sizes and shapes can be generated. Although blocks of 
hypothetical layers 404 are depicted in Fig. 4 as including three hypothetical layers 402, 
it should be recognized that they can include any number of hypothetical layers 402. 

[0028] Diffraction calculations, which characterize the behavior of diffraction 

signals, are generated for each hypothetical layer 402 within a block of hypothetical 
layers 404. Diffraction calculations for each block of hypothetical layers 404 are then 
generated by aggregating the diffraction calculations for each hypothetical layer 402 
within each block of hypothetical layers 404. The diffraction calculations and the blocks 
of hypothetical layers 404 are then stored. More particularly, in one exemplary 
embodiment, pairs of diffraction calculations and blocks of hypothetical layers 404 are 
stored in a cache 406. 

[0029] For example, for block of hypothetical layers 404.1, diffraction 

calculations are generated for hypothetical layers 402.1, 402.2, and 402.3. A diffraction 
calculation for block of hypothetical layers 404.1 is then generated by aggregating 
diffraction calculations for hypothetical layers 402.1, 402.2, and 402.3. Block of 
hypothetical layers 404.1 and the diffraction calculations associated with block of 
hypothetical layers 404.1 are then stored in cache 406. In a similar manner, diffraction 



calculations for blocks of hypothetical layers 404.2, 404.3, and 404.4 are generated and 
stored in cache 406. 



[0030] Although four blocks of hypothetical layers 404 are depicted and 

described as being generated and stored in cache 406, it should be recognized that any 
number of blocks of hypothetical layers 404 of various shapes and configurations can be 
generated and stored in cache 406. Indeed, in use, cache 406 may contain tens and 
hundreds of thousands of blocks of hypothetical layers 404 and diffraction calculations. 

[0031] Additionally, although Fig. 4 depicts blocks of hypothetical layers 404 

being stored in cache 406 in a graphical format, blocks of hypothetical layers 404 can be 
stored in various formats. For example, blocks of hypothetical layers 404 can be stored 
using parameters that define their shapes. 

[0032] In one exemplary embodiment, simulated-diffraction signals for a set of 

hypothetical profiles to be stored in library 185 can be generated based on the diffraction 
calculations for blocks of hypothetical layers 404 stored in cache 406. More particularly, 
for each hypothetical profile to be stored in library 185, one or more blocks of 
hypothetical layers 404 that characterize the hypothetical profile are selected from those 
stored in cache 406. Boundary conditions are then applied to generate the simulated- 
diffraction signal for the hypothetical profile. The simulated-diffraction signal and the 
hypothetical profile are then stored in library 185. The hypothetical profile can be stored 
in various formats, such as graphically, using the parameters that define the hypothetical 
profile, or both. 

[0033] In some applications, one block of hypothetical layers 404 can be used to 

characterize a hypothetical profile. For example, with reference to Fig. 5A, assume for 
the sake of this example that a simulated-diffraction signal is to be generated for 
hypothetical profile 3 00 A. As described above, a block of hypothetical layers 404 from 
cache 406 is selected that characterizes hypothetical profile 300A. The appropriate block 
of hypothetical layers 404 can be selected using an error minimization algorithm, such as 
a sum-of-squares algorithm. 
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[0034] As depicted in Fig. 5 A, in this example, block of hypothetical layers 404. 1 

is selected from cache 406. The diffraction calculation associated with block of 
hypothetical layers 404.1 is then retrieved from cache 406. Boundary conditions are then 
applied to generate the simulated-diffraction signal for hypothetical profile 300A. The 
simulated-diffraction signal and hypothetical profile 300A are then stored in library 185. 

[0035] As described above, the parameters that define hypothetical profile 300A 

can be varied to define another hypothetical profile. In this example, with reference to 
Fig. 5B, assume that the parameters are varied to define hypothetical profile 300B. 

[0036] As depicted in Fig. 5B, block of hypothetical layers 404.4 is selected from 

cache 406. The diffraction calculation associated with block of hypothetical layers 404.4 
is then retrieved from cache 406. Boundary conditions are then applied to generate the 
simulated-diffraction signal for hypothetical profile 300B. The simulated-diffraction 
signal and hypothetical profile 300B are then stored in library 185. 

[0037] In this manner, any number of hypothetical profiles and simulated- 

diffraction signals can be generated and stored in library 185 using blocks of hypothetical 
layers 404 and the diffraction calculations in cache 406. As described above, in the 
present exemplary embodiment, blocks of hypothetical layers 404 and the diffraction 
calculations in cache 406 are generated before generating library 185. As such, the 
calculations needed to generate the diffraction calculations associated with blocks of 
hypothetical layers 404 do not need to be performed at the time that library 185 is 
generated. This can have the advantage that the amount of time and calculations needed 
to generate library 185 can be reduced. 

[0038] In some applications, multiple blocks of hypothetical layers 404 can be 

used to characterize a hypothetical profile. For example, with reference to Fig. 6A, 
profile 300A can be characterized by blocks of hypothetical layers 404.1, 404.2, and 
404.3. In the present embodiment, an error minimization algorithm can be utilized to 
determine the appropriate blocks of hypothetical layers 404 to use in characterizing 
profile 300A. 
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[0039] In the present example, for profile 300A, the diffraction calculations 

associated with blocks of hypothetical layers 404. 1, 404.2, and 404.3 are retrieved from 
cache 406. Boundary conditions are then applied to generate the simulated-diffraction 
signal for profile 300A. More particularly, the boundary conditions at the top and bottom 
of profile 3 00 A are applied. The resulting simulated-diffraction signal and profile 300A 
are then stored in library 185. 

[0040] As described above, the parameters that define hypothetical profile 300A 

can be varied to define another profile. In this example, with reference to Fig. 6B, 
assume that the parameters are varied to define hypothetical profile 300B. 

[0041 ] As depicted in Fig. 6B, hypothetical profile 300B can be characterized by 

blocks of hypothetical layers 404.1, 404.3, and 404.4. As such, the diffraction 
calculations associated with blocks of hypothetical layers 404.1, 404.3, and 404.4 are 
retrieved from cache 406. Boundary conditions are then applied to generate the 
simulated-diffraction signal for hypothetical profile 300B. The resulting simulated- 
diffraction signal and hypothetical profile 300B are then stored in library 185. 

[0042] In some applications, a hypothetical profile can include multiple materials. 

For example, with reference to Fig. 7 A, assume that hypothetical profile 300 is used to 
characterize an actual profile of a grating periodic formed from two materials. More 
particularly, as depicted in Fig. 7A, hypothetical profile 300 includes a first layer 702 and 
a second layer 704. In the present example, assume that first layer 702 represents an 
oxide layer of the actual profile, and second layer 704 represents a metal layer of the 
actual profile. Note that hypothetical profile 300 can include any number of layers to 
represent layers of any number of materials in an actual profile. 

[0043] With reference to Fig. 7B, in the present example, assume that cache 406 

includes blocks of hypothetical layers 404 of various materials as well as various shapes. 
For example, assume that blocks of hypothetical profiles 404.1 and 404.3 represent oxide 
layers, and hypothetical profiles 404.2 and 404.4 represent metal layers. As such, as 
depicted in Fig. 7B, profile 300 can be characterized by blocks of hypothetical layers 
404.3 and 404.4 from cache 406. Boundary conditions are then applied to generate the 



12 



simulated-diffraction signal for profile 300. The simulated-diffraction signal and profile 
300 are then stored in library 185. 

[0044] Thus far, hypothetical layers 402 and blocks of hypothetical layers 404 

have been depicted as having rectangular and trapezoidal shapes, respectively. As will be 
described in greater detail below, the diffraction calculation for a hypothetical layer 402 
depends on its width but not on its height. However, the diffraction calculation for a 
block of hypothetical layers 404 depends on its height as well as its width. Therefore, the 
blocks of hypothetical layers 404 stored in cache 406 are characterized and indexed, in 
part, by their width and their height. More particularly, when blocks of hypothetical 
layers 404 have symmetric-trapezoidal shapes, they can be indexed by their height, 
bottom width (bottom CD), and top width (top CD). However, as noted earlier, blocks of 
hypothetical layers 404 can have various shapes. As such, they can be characterized and 
indexed using any number of parameters. 

[0045] With reference to Fig. 8, in one exemplary embodiment, library 185 can 

be generated using computer system 800. As noted above and as will be described in 
greater detail below, the process of generating library 185 can involve performing a large 
number of complex calculations. 

[0046] For example, assume that a set of hypothetical profiles is to be generated 

for a periodic grating that is to have a profile with a bottom CD of 200 nm. Assume that 
a 10% process variation is expected for the bottom CD of the periodic grating, which 
means that the bottom CD is expected to vary between about 180 to about 220 nm. Also 
assume that the top CD is expected to vary between about 160 to about 180 nm. Assume 
that the nominal thickness (i.e., the height) of the periodic grating is to be about 500 nm, 
and that a 10% process variation is expected, which means that the height can vary 
between about 450 to about 550 nm. Now assume that the desired resolution is 1 nm, 
which means that each parameter of the hypothetical profiles is varied by an increment of 
1 nm. 

[0047] In generating the set of hypothetical profiles, the top CD of the 

hypothetical profiles is varied between 160 to 180 nm in steps of 1 nm. The bottom CD 
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of the hypothetical profiles is varied between 180 to 220 nm in steps of 1 nm. The 
thickness/height of the hypothetical profiles is varied between 450 to 550 nm in steps of 1 
nm. Thus, in this example, a total of 87,000 hypothetical profiles (i.e., 21 variations of 
the top CD multiplied by 41 variations of the bottom CD multiplied by 101 variations of 
the thickness/height). 

[0048] Now, assume that diffraction calculations are to be generated for each 
hypothetical profile at 53 different wavelengths. Also assume that each diffraction 
calculation uses 6 matrices with 9 orders and 8 bytes, which totals 17 kbytes. As such, in 
this example, to store the diffraction calculations for all of the 87,000 hypothetical 
profiles at 53 different wavelengths, a total of 78 Gigabytes is needed. 

[0049] As such, computer system 800 can include multiple processors 802 

configured to perform portions of the computations in parallel. However, computer 
system 800 can be configured with a single processor 802. 

[0050] Additionally, computer system 800 can include a memory 804 configured 

with a large amount of memory, such as 8 Gigabytes, 16 Gigabytes, 32 Gigabytes, and 
the like, that can be accessed by the multiple processors 802. It should be recognized, 
however, computer system 800 can be configured with any number and size of memories 
804. 

[0051] In one exemplary embodiment, as depicted in Fig. 8, cache 406 can reside 

in memory 804. As such, the blocks of hypothetical layers and hypothetical profiles 
stored in cache 406 can be more quickly accessed by computer system 800, and more 
particularly processors 802, than if cache 406 resided on a hard drive. 

[0052] Additionally, in one exemplary embodiment, library 1 85 can reside on 
various computer-readable storage media. For example, library 185 can reside on a 
compact disk that is written to by computer system 800 when library 1 85 is generated, 
and read by signal processing module 190 (Fig. 1) when library 185 is used to determine 
the profile of periodic grating 145 (Fig. 1) on wafer 140 (Fig. 1). 
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3. Rigorous Coupled- Wave Analysis for Incident Radiation 

[0053] As described above, simulated-diffraction signals are generated for 

hypothetical profiles based on diffraction calculation generated for hypothetical layers 
402 and blocks of hypothetical layers 404. As will be described below, in one exemplary 
embodiment, simulated-diffraction signals and diffraction calculations can be generated 
by applying Maxwell's equations and using a numerical analysis technique to solve 
Maxwell's equations. More particularly, in the exemplary embodiment described below, 
rigorous coupled-wave analysis is used. It should be noted, however, that various 
numerical analysis techniques, including variations of the rigorous coupled-wave analysis 
technique described below, can be used. 

[0054] With reference to Fig. 9, a section of periodic grating 145 is shown. In 

Fig. 9, the section of the grating 145 is depicted as having three ridges 902, which have a 
triangular cross-section. It should be noted that ridges can have shapes that are 
considerably more complex, and the categories of "ridges" and "troughs" may be ill- 
defined. For the sake of convenience and clarity, the term "ridge" will be used hereafter 
to refer to one period of a periodic structure on a substrate. Each ridge 902 is considered 
to extend infinitely in the +y and -y directions, and an infinite, regularly-spaced series of 
such ridges 902 are considered to extend in the +x and -x directions. The ridges 902 are 
atop a deposited film 904, and the film 904 is atop a substrate 906, which is considered to 
extend semi-infinitely in the +z direction. The normal vector n to the grating is in the -z 
direction. 

[0055] Fig. 9 illustrates the variables associated with a mathematical analysis of a 

diffraction grating according to the present exemplary embodiment. In particular: 

• 9 is the angle between the Poynting vector 908 of the 
incident electromagnetic radiation 910 and the normal vector nof the grating 145. 
The Poynting vector 908 and the normal vector n define the plane of incidence 
912. 

• <|> is the azimuthal angle of the incident electromagnetic 
radiation 910, i.e., the angle between the direction of periodicity of the grating, 
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which in Fig. 9 is along the x axis, and the plane of incidence 912. For ease of 
presentation, in the mathematical analysis of the present specification the 
azimuthal angle <|> is set to zero. 

• i)/ is the angle between the electric-field vector E of the 
incident electromagnetic radiation 910 and the plane of incidence 912, i.e., 
between the electric field vector E and its projection E } on the plane of incidence 
912. When (j)=0 and the incident electromagnetic radiation 910 is polarized so 
that V|/=7i/2 5 the electric-field vector £is perpendicular to the plane of incidence 
912 and the magnetic-field vector H lies in the plane of incidence 912, and this is 
referred to as the TE polarization. When <[)=0 and the incident electromagnetic 
radiation 910 is polarized so that \|/=0, the magnetic-field vector His 

perpendicular to the plane of incidence 912 and the electric-field vector i?lies in 
the plane of incidence 912, and this is referred to as the TM polarization. Any 
planar polarization is a combination of in-phase TE and TM polarizations. The 
method of the present invention described below can be applied to any 
polarization that is a superposition of TE and TM polarizations by computing the 
diffraction of the TE and TM components separately and summing them. 

Furthermore, although the 'off-axis' ty£0 case is more complex because it cannot 
be separated into TE and TM components, the present invention is applicable to 
off-axis incidence radiation as well. 

• X is the wavelength of the incident electromagnetic 

radiation 910. 

[0056] Fig. 10A shows a cross-sectional view of two ridges 902 of periodic 

grating 145 illustrating the variables associated with a mathematical description of the 
dimensions of the diffraction grating 145. In particular: 

• L is the number of the layers into which the system is 
divided. Layers 0 and L are considered to be semi-infinite layers. Layer 0 is an 
"atmospheric" layer 1002, such as vacuum or air, which typically has a refractive 
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index no of unity. Layer L is a "substrate" layer 906, which is typically silicon or 
germanium in semiconductor applications. In the case of the exemplary grating 
145 of Fig. 10A, the grating 145 has ten layers with the atmospheric layer 1002 
being the zeroeth layer 1025.0, the ridges 902 being in the first through seventh 
layers 1025.1 through 1025.7, the thin film 904 being the eighth layer 1025.8, and 
the substrate 906 being the ninth layer 1025.9. For the mathematical analysis 
described below, the thin-film 904 is considered as a periodic portion of the ridge 
902 with a width d equal to the pitch D. The portion of ridge 902 within each 
intermediate layer 1025.1 through 1025.(L-1) is approximated by a thin planar 
slab 1026 having a rectangular cross-section. Genetically or collectively, the 
layers are assigned reference numeral 1025, and, depending on context, "layers 
1025" may be considered to include the atmospheric layer 1002 and/or the 
substrate 906. In general, any geometry of ridges 902 with a cross-section which 
does not consist solely of vertical and horizontal sections can be better 
approximated using a large number of layers 1025. 

• D is the periodicity length or pitch, i.e., the length between 
equivalent points on pairs of adjacent ridges 902. 

• di is the width of the rectangular ridge slab 1026./ in the /th 

layer 1025./. 

• U is the thickness of the rectangular ridge slab 1026./ in the 
/th layer 1025./ for 1</<(Z-1). The thicknesses t\ of the layers 1025 are chosen 
such that every vertical line segment within a layer 1025 passes through only a 
single material. For instance, if in Fig. 10 the materials in layers 1025.4, 1025.5, 
and 1025.6 are the same, but different than the materials in layers 1025.3 and 
1025.7, than it would be acceptable to combine layers 1025.4 and 1025.5, or 
layers 1025.5 and 1025.6, or layers 1025.4, 1025.5 and 1025.6 into a single layer. 
However, it would not be acceptable to combine layers 1025.3 and 1025.4, or 
layers 1025.6 and 1025.7 into a single layer. 

• ni is the index of refraction of the material in the 
rectangular ridge slab 1026 of the /th layer 1025./. 
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[0057] In determining the diffraction generated by grating 145, a Fourier space 

version of Maxwell's equations is used. As shown in the calculation process flow 
diagram of Fig. 1 1, the permittivities z\(x) for each layer / are determined or acquired 
(1 1 10), and a one-dimensional Fourier transformation of the permittivity Zj(x) of each 
layer / is performed (1 1 12) along the direction of periodicity, x , of the periodic grating 
145 to provide the harmonic components of the permittivity Zi th where i is the order of the 
harmonic component. In Figs. 1 1, 12, 13, and 14 process steps are shown enclosed 
within ovals or rectangles with rounded corners, and the results of calculations are shown 
enclosed within rectangles with sharp corners. When appropriate in Fig. 11, equation 
numbers are used in lieu of, or in addition to, reference numerals. In particular, the real- 
space permittivity z\(x) of the /th layer is related to the permittivity harmonics sij of the 
Ith layer by 

j=-O0 

[0058] Therefore, via the inverse transform, 

s i)0 =^f + «od-f), (1-1.2) 

and for i not equal to zero, 

sm(m — ) 

s u =(K-nl) P- e -J*» t (1.1.3) 

711 

where n r is the index of refraction of the material in the ridges 902 in layer /, the index of 
refraction n 0 of the atmospheric layer 1002 is typically near unity, and p is the x-offset of 
the center of the central rectangular ridge slab 1026./ (i.e., the ridge 902 nearest x=0 9 
where generally it is attempted to position the x=0 point at the center of a ridge 902) from 
the origin. The present specification explicitly addresses periodic gratings where a single 
ridge material and a single atmospheric material are found along any line in the x- 
direction. However, gratings can have more than one ridge material along a line in the x- 
direction. 

[0059] In the present embodiment, it is convenient to define the (2o+l) x (2o+l) 

Toeplitz-form, permittivity harmonics matrix E\ as 
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°/,0 


Si , 


°l-2 


°/,-2o 


*/.! 






^/ ,-(2o-l) 








£ l,-(2o-2) 


£ l,2o 


S l,{2o-X) 


S l,{2o-2) 


S /,0 



(1.1.4) 



[0060] As will be seen below, to perform a TE-polarization calculation where 

oth-order harmonic components of the electric field E and magnetic field H are used, it 
is necessary to use harmonics of the permittivity s u up to order 2o. 



[0061] For the TE polarization, in the atmospheric layer the electric field E is 

formulated (1124) as 

E Q y = exp(-jk 0 n 0 (sin Be + cosBz) + J] R. vq>(-j(k xi x - k 0zi z)) , (1 .2. 1) 

where the term on the left of the right-hand side of equation (1 .2.1) is an incoming plane 
wave at an angle of incidence 6, the term on the right of the right-hand side of equation 
(1 .2.1) is a sum of reflected plane waves and i? z is the magnitude of the ith component of 
the reflected wave, and the wave vectors ko and (k xi , ko tZ d are given by 



K = ^ = <D(/i 0 e 0 ) 1/2 , 
X 

k » =K(n o sm(0) -K—)) 



(1.2.2) 
(1.2.3) 



and 



K, 2i = ■ 



(1.2.4) 



[-jk 0 (n?-(kJK) 2 y' 2 
where the value of ko fZ i is chosen from equation (1 .2.4), i.e., from the top or the bottom of 
the expression, to provide Re(^ z /) - Im(£^) > 0. This insures that k\ zi has a positive 

real part, so that energy is conserved. It is easily confirmed that in the atmospheric layer 
1002, the reflected wave vector (k xi , ko, z d has a magnitude equal to that of the incoming 
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wave vector k 0 n 0 . The magnetic field H in the atmospheric layer 1002 is generated from 
the electric field Eby Maxwell's equation (1.3.1) provided below. 

[0062] The x-components k xi of the outgoing wave vectors satisfy the Floquet 

condition (which is also called Bloch's Theorem, see Solid State Physics, N.W. Ashcroft 
and N.D. Mermin, Saunders College, Philadelphia, 1976, pages 133-134) in each of the 
layers 1025 containing the periodic ridges 902, and therefore, due to the boundary 
conditions, in the atmospheric layer 1002 and the substrate layer 906 as well. That is, for 
a system having an ^-dimensional periodicity given by 

/(0 = /(r>I^) 5 (1.2.5) 

where d. are the basis vectors of the periodic system, and m z takes on positive and 
negative integer values, the Floquet condition requires that the wave vectors k satisfy 

_ _ J* _^ 

£ = £ 0 +27i£ wl A'> (1.2.6) 
1=1 

where b { are the reciprocal lattice vectors given by 

{b r dj)=* t , (1.2.7) 
k 0 is the wave vector of a free-space solution, and 5~ is the Kronecker delta function. In 
the case of the layers 1025 of the periodic grating 145 of Figs. 9 and 10, which have the 
single reciprocal lattice vector b is x I D , thereby providing the relationship of equation 
(1.2.3). 

[0063] It may be noted that the formulation given above for the electric field in 

the atmospheric layer 1002, although it is an expansion in terms of plane waves, is not 
determined via a Fourier transform of a real-space formulation. Rather, the formulation 
is produced (1 124) a priori based on the Floquet condition and the requirements that both 
the incoming and outgoing radiation have wave vectors of magnitude n^. Similarly, the 
plane wave expansion for the electric field in the substrate layer 906 is produced (1 124) a 

priori. In the substrate layer 906, the electric field E is formulated (1 124) as a 
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transmitted wave which is a sum of plane waves where the x-components k xi of the wave 
vectors (k xif k 0tZ i) satisfy the Floquet condition, i.e., 



K =y L T t «p(-i(*^ + KA Z - > 0-2-8) 

i /-I 

where 



*o(^-(* rf /*o) 2 ) l/2 

{-jk Q {n 2 L -{kjk 0 yy 



where the value ofk LtZi is chosen from equation (L2.9), i.e., from the top or the bottom of 
the expression, to provide Re(k LfZi ) - lm(k LtZi ) > 0, insuring that energy is conserved. 

[0064] The plane wave expansions for the electric and magnetic fields in the 

intermediate layers 1025.1 through 1025.(Z-1) are also produced (1 134) a priori based on 

the Floquet condition. The electric field E l y in the /th layer is formulated (1 134) as a 

plane wave expansion along the direction of periodicity, x , i.e., 

i 

where S lyi (z) is the z-dependent electric field harmonic amplitude for the /th layer and 

the zth harmonic. Similarly, the magnetic field H l y in the /th layer is formulated (1 134) 
as a plane wave expansion along the direction of periodicity, x , i.e., 

H l>x = -7X^) 1/2 S^(^)exp(- 7 ^) , (1.2.11) 
Ho / 

where U lxi (z) is the z-dependent magnetic field harmonic amplitude for the /th layer and 
the zth harmonic. 

[0065] According to Maxwell's equations, the electric and magnetic fields within 

a layer are related by 
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#,=(— )Vx£, , (1.3.1) 
<yu 0 

and 

E,=( ~ J WxH, . (1.3.2) 

Applying (1142) the first Maxwell's equation (1.3.1) to equations (1.2.10) and (1.2.11) 
provides a first relationship between the electric and magnetic field harmonic amplitudes 
Si and Ui of the /th layer: 



gMf> 



Similarly, applying (1 141) the second Maxwell's equation (1.3.2) to equations (1.2.10) 
and (1.2.1 1), and taking advantage of the relationship 

2nh _ 

^ja + = *x(/-A) (1 .3.4) 

which follows from equation (1.2.3), provides a second relationship between the electric 
and magnetic field harmonic amplitudes Si and U\ for the /th layer: 



dU, 

l,x 



^Uyi X S (i-p) ^Uyp ' (1.3.5) 



While equation (1.3.3) only couples harmonic amplitudes of the same order i, equation 
(1.3.5) couples harmonic amplitudes Si and Ui between harmonic orders. In equation 
(1.3.5), permittivity harmonics e,- from order -2o to +2o are required to couple harmonic 
amplitudes Si and Ui of orders between -o and +o. 

[0066] Combining equations (1 .3.3) and (1 .3.5) and truncating the calculation to 

order o in the harmonic amplitude S provides (1 145) a second-order differential matrix 
equation having the form of a wave equation, Le., 
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= Mk y ] , (1-3.6) 

z ' = k 0 z, the wave- vector matrix [A{\ is defined as 

[4] = feJ-fe] , (1-3-7) 

where [K x ] is a diagonal matrix with the element being equal to (k xi /ko), the 
permittivity harmonics matrix [Ei\ is defined above in equation (1.1 .4), and [ Si, y ] and 
are column vectors with indices i running from ~o to +o, i.e., 

ky} = 

[0067] By writing (1 1 50) the homogeneous solution of equation (1 .3 .6) as an 

expansion in pairs of exponentials, i.e., 



d 2 S l>y 

dz' 2 

_ 



Sl,y,o 



(1.3.8) 



S l,yi ( Z ) = Z Vhm eXp(-*0?/,«^) + C2 l,m ^MK^m (* " h ))] - 0 -3.9) 

l«=l 

its functional form is maintained upon second-order differentiation by z thereby taking 
the form of an eigenequation. Solution (1 147) of the eigenequation 

UMh bM]> (i-3.io) 

provides (1 148) a diagonal eigenvalue matrix [x/] formed from the eigenvalues x/ jYn of the 
wave- vector matrix [Aj\, and an eigenvector matrix [Wi] of entries Wij tm , where Wij ftn is the 
zth entry of the mth eigenvector of [A{\. A diagonal root-eigenvalue matrix [Q{] is defined 
to be diagonal entries qu which are the positive real portion of the square roots of the 
eigenvalues x/ The coefficients cl and c2 are, as yet, undetermined. 

[0068] By applying equation (1 .3.3) to equation (1 .3.9) it is found that 
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2o+l 



Ut, xi (z) = Z v u« I" ex P(-^o?/^^) + ^ exp(£ 0 ^ m (z - *,))] (1.3.1 1) 



where = ^ m w/^; w . The matrix [ Vj], to be used below, is composed of entries v/ >z >. 

[0069] The coefficients cl and c2 in the homogeneous solutions of equations 

(1 .3.9) and (1.3.1 1) are determined by applying (1 155) the requirement that the tangential 
electric and magnetic fields be continuous at the boundary between each pair of adjacent 
layers 1025.//(1025.(M-1). At the boundary between the atmospheric layer 1002 and the 
first layer 1025. 1, continuity of the electric field E y and the magnetic field H x requires 



jn 0 cos(0)5 /o 



/ 



R = 



w x w x x x 
Y\ -v x x x 



cl x 



(1.4.1) 



where Y 0 is a diagonal matrix with entries (ko rZ i/k 0 ), Xi is a diagonal layer-translation 
matrix with elements cxp(-ko qi m ti), R is a vector consisting of entries from R. 0 to R +0 
and c7/and c2i are vectors consisting of entries from cli t o to cli t 2o+n and c2 it0 to c2i >2o +n 
respectively. The top half of matrix equation (1.4.1) provides matching of the electric 
field E y across the boundary of the atmospheric layer 1025.0 and the first layer 1025.1, 
the bottom half of matrix equation (1.4.1) provides matching of the magnetic field H x 
across the layer boundary 1025.0/1025.1, the vector on the far left is the contribution 
from the incoming radiation 910 in the atmospheric layer 1002, the second vector on the 
left is the contribution from the reflected radiation 914 in the atmospheric layer 1002, and 
the portion on the right represents the fields E y and//* in the first layer 1025.1. 

[0070] At the boundary between adjacent intermediate layers 1 025./ and 

1025.(/+1), continuity of the electric field E y and the magnetic field H x requires 











W, 






7m*m 




_c2 M _ 




J, 


-v,x,_ 


[c2,_ 



(1.4.2) 



where the top and bottom halves of the vector equation provide matching of the electric 
field E y and the magnetic field H x , respectively, across the I A II layer boundary. 
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[0071] At the boundary between the (I-l)th layer 1025.(L-1) and the substrate 
layer 906, continuity of the electric field E y and the magnetic field H x requires 





















J Y L. 



where, as above, the top and bottom halves of the vector equation provides matching of 
the electric field E y and the magnetic field H x , respectively. In contrast with equation 
(1 Al), there is only a single term on the right since there is no incident radiation in the 
substrate 906. 

[0072] The method (1200) of calculation for the diffracted reflectivity of TM- 

polarized incident electromagnetic radiation 910 shown in Fig. 12 parallels the method 
(1 100) described above and shown in Fig. 1 1 for the diffracted reflectivity of TE- 
polarized incident electromagnetic radiation 910. The variables describing the geometry 
of the grating 145 and the geometry of the incident radiation 910 are as shown in Figs. 9 
and 10A. However, for TM-polarized incident radiation 910 the electric field vector 2? is 
in the plane of incidence 912, and the magnetic field vector H is perpendicular to the 
plane of incidence 912. The similarity in the TE- and TM-polarization RCWA 
calculations, in part, motivates the use of the term "electromagnetic field" in the present 
specification to refer generically to either or both the electric field and/or the magnetic 
field of the electromagnetic radiation. 

[0073] The most substantial difference between the calculation (1200) for TE- 

polarized incident electromagnetic radiation 910 and the calculation (1 100) for TM- 
polarized incident electromagnetic radiation 910 is that for TM-polarized incident 
radiation the wave-vector matrix [A{\ may be defined according to any of three 
formulations: 

M=[E t ]tK,][P,][K x ]-[rj), (1.5.1) 

or 

[A,]= [PJ- 1 ([KJiE^iKJ-il]) , (1.5.2) 
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or 



[4]=[iy([^j[£j-^j-m)- (i.5.3) 

In the above equations, the inverse-permittivity harmonics matrix Pi is defined as 



^1,0 




*/.-2 • 


" n i,-2o 




71 1,0 




K l,-(2o-\) 




71 1,1 




n i-(2o-2) 


*l.2o 




^1,(20-2) ' 


1*1,0 



(1.5.4) 



where 



1 ^ ,.2nh s 
e,(x) h TL D 



(1.5.5) 



so via the inverse Fourier transform, 



1 d i I/, d i^ 
n„ JJ n 0 D 



(1.5.6) 



and for h not equal to zero, 



d, 

i i sin^-A 
1 - D 0 -mp/d 



n 2 r n\ ' nh 



(1.5.7) 



The fourth possible formulation for the wave- vector matrix [A{\ 

W=[P l T 1 ([K x ][P l ][K x ]-[I]) , 



(1.5.8) 



tends not to be computationally useful. Based on these definitions of the wave- vector 
matrix [A{\, the matrix [Vj\ is defined as 



[Vi\ = [E,]- 1 [W][Qj] 
when [A] is defined as in equation (1.5.1), 

[Vi\ = [Pi\[Wi\[Q{\ 
when [A] is defined as in equation (1.5.2), and 



(1.5.9) 



(1.5.10) 
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[Vii^[Ei\ A [Wi\[Qi\ (1.5.11) 

when [^4] is defined as in equation (1.5.3). The convergence of the RCWA calculation 
will vary depending on the formulation for the wave- vector matrix [A{\ 9 and the 
particulars of the grating 145 and the incident radiation 910. 

[0074] Following through the calculation (1200) for TM-polarized incident 

electromagnetic radiation 910 paralleling the calculation (1 100) described above for TE- 
polarized incident electromagnetic radiation 910, it is found that the reflectivity vector R 
is related to the coefficients cl \ and c2\ by 



i0 



jcos(0)5 iO /n o 
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R = 




w x x x ~ 






_-J Z o_ 




7x 


-v x x x _ 





(1.6.1) 



where Z 0 is a diagonal matrix with entries (ko, zi /n 0 2 k 0 \ and, as above, X\ is the diagonal 
layer-translation matrix with elements exp(-£ 0 q Lm tj). Similarly, following through the 
calculation (1200) for TM-polarized incident electromagnetic radiation 910 paralleling 
the calculation (1 100) described above for TE-polarized incident electromagnetic 
radiation 910, it is found that the coefficients cl/ and cl\ are related to the coefficients cl/_ 
i and c2\-i by 







"cl H " 




>/ 


W t X,' 








_c2 M _ 




J i 


-v t x tm 


ft] 



(1.6.2) 



and the transmission vector J is related to the coefficients c\ L .\ and cl L .\ by 







' I ' 






J Z L. 



(1.6.3) 



[0075] Equation (1 .6.1) only differs from equation (1 .4.1) by the substitution of 

Z 0 for Y 0 and the replacement of no with l/n 0 on the left side of the equation. Similarly, 
equation (1 .6.3) only differs from equation (1 .4.3) by the substitution of Z 0 for Y 0 on the 
right side of the equation. It is important to note that equation (1 .6.2) is the same as 
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equation (1.4.2), the right side of equation (1.6.1) is the same as the right side of equation 
(1.4.1), and the left side of equation (1.6.3) is the same as the left side of equation (1.4.3). 

[0076] As disclosed in pending patent application serial number 09/770,997, filed 

January 25, 2001, entitled Caching of Intra-Layer Calculations for Rapid Rigorous 
Coupled-Wave Analyses, by the same inventors, and incorporated herein by reference, 
for TE-polarized incident electromagnetic radiation 910, matrix equation (1.4.1), matrix 
equation (1.4.3), and the (L-l) matrix equations (1.4.2) can be combined (1 160) to 
provide the system matrix equation: 
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w x x x 
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-w L . x x L _ x 
-v L .,x L _ x 


-w L . x 
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T 




0 



and this boundary-matched system matrix equation may be solved to provide the 
reflectivity R t for each harmonic order i. For TE-polarized incident electromagnetic 
radiation 910, equations (1.6.1), (1.6.2), and (1.6.3) maybe solved via a similar system 
matrix equation. 

4. Transfer Matrix Formulation of the System Equation 

[0077] As taught by "Stable Implementation of the Rigorous Coupled- Wave 

Analysis for Surface-Relief Dielectric Gratings: Enhanced Transmittance Matrix 
Approach", E. B. Grann and D. A. Pommet, Opt Soc. Am. A, vol. 12, 1077-1086, May 
1995, equation (1.4.2) may be written as 



"cl M " 






-1 




w x x x - 












J, 


-v,x,_ 


c2 l 



and, for TM-polarized incident electromagnetic radiation 910, the reflectivity R may be 
coupled to the transmission Tby repeated substitution of equation (2.1.1) for values of / 
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from 1 to (1-1) into equations (1.4.1) and (1.4.3) to provide (1 160) system matrix 
equation 



ljn 0 cos(0)5 i0 }l-jY 0 \ R 0 < 



L-\ 



y -r,x t lr,x, -v t 



' i 



(2.1.2) 

where the pairs of matrices in curly brackets are arranged in ascending order in the layer 
variable / from left to right. Equation (2.1.2) may be solved to provide the reflectivity R. 
Similarly, for TE-polarized incident electromagnetic radiation 910 the reflectivity R may 
be coupled to the transmission T to provide (1 160) system matrix equation 



T , 
(2.1.3) 
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which may be solved for the reflectivity R. 

[0078] According to the present specification, the matrix 





w t x,- 






7i 


-v,x,_ 


y l x l 





(2.1.4) 



is termed the "transfer matrix" for the / th layer, and equations (2. 1 .2) and (2. 1 .3) may be 
written as 
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> 0 cos(^)5. 0 _ 




.-J Y o. 
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J Y L. 
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(2.1.5) 



and 
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Jcos(0)S io /n o 
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J Z L. 


T , 



(2.1.6) 

respectively. The transfer matrix [Ti\ is termed the 'transmittance' matrix in some 
references, such as the above-cited reference by E.B. Grann and D.A. Pommet. Note that 
in the present exemplary embodiment just as the wave-vector matrix [A{\, and its 
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eigenvector matrix [Wi\ and root-eigenvalue matrix [Qi] are only dependent on intra-layer 
parameters for the / th layer (and incident-radiation parameters), and the layer-transfer 
matrix Ti is also only dependent on intra-layer parameters for the 1th layer (and incident- 
radiation parameters). However, in contrast with the wave-vector matrix \A{\ 9 and its 
eigenvector matrix [Wi] and root-eigenvalue matrix [Qj\ 9 the layer-translation matrix [Xi] 
which appears in the layer-transfer matrix [T/] is dependent on the thickness */ of the / th 
layer. Although thickness t is an intra-layer parameter, for reasons that will become 
apparent below, thickness t will be classified as an aggregate-layer parameter in the 
present specification. 

[0079] In the present exemplary embodiment, the "aggregate-layer transfer 

matrix" for layers li to Ij can be defined as: 



l=li 



v, 



\ w l x l Tw l x l w x 



(2.1.7) 



where the pairs of matrices on the right side of equation (2.1.7) are arranged in ascending 
order in the layer variable / from left to right. For example, 



[r^ 4 ] = 



>3 


w 3 x 3 ~ 


~w 3 x 3 


w 3 ~ 


-1 


>4 


w 4 x 4 ~ 


~w 4 x 4 


w 4 ~ 


J 3 


-V 3 X 3 _ 


V 3 X 3 








-v 4 x 4 _ 


y<x 4 





(2.1.8) 



Note that the aggregate-layer transfer matrix [rV»/ 7 ] is only dependent on the 
characteristics of the /z th through if 1 layers (and the incident-radiation parameters). In 
terms of the aggregate-layer transfer matrices [T/i->y], the system matrix equations for a 
profile which is discretized into n+l aggregate-layer trapezoids are 



jn 0 cos(6>)5 
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(2.1.9) 
T , 

(2.1.10) 
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for TM- and TE-polarizations, respectively. As described in detail below, caching and 
retrieval of aggregate-layer transfer matrices [r^j] can be used to reduce the 
computation time required for generating spectra using rigorous coupled-wave analyses. 

5. Caching and Retrieval of Aggregate-Layer Transfer Matrices 

[0080] As presented above, the calculation of the diffraction of incident TE- 

polarized or TM-polarized radiation 910 from a periodic grating involves the generation 
of a system matrix equation (2.1.5) or (2.1.6), respectively, and its solution. The most 
computationally expensive portion of the RCWA computations flowcharted in Figs. 1 1 
and 12 can be the solution (1 147) and (1247) for the eigenvectors and eigenvalues Ty 
of wave-vector matrix [Ai] of equation (1.3.7). The accuracy of the calculation of the 
eigenvectors wij and eigenvalues is dependent on the number of orders o utilized. As 
the number of orders o is increased, the computation time for solving the eigensystem 
increases exponentially with o. When performed in a typical computing environment 
with o=9 harmonic orders, the calculation of the eigenvectors and eigenvalues can take 
more than 85% of the total computation time. 

[0081 ] In the present exemplary embodiment, a computer system 1 500 can be 

used to perform the various calculations and operations described above. As depicted in 
Fig. 15, computer system 1500 can include an information input/output (I/O) equipment 
1505, which is interfaced to a computer 1510 which includes a central processing unit 
(CPU) 1515 and a memory 1520. Although the CPU 1515 is depicted in Fig. 15 as a 
single unit, the CPU 1515 can include a number of processors configured to process in 
parallel. In order of increasing preference, the CPU 1515 includes four, eight, sixteen, 
thirty-two, or more high-speed parallel processors. The parallel processing capabilities of 
the CPU 1515 of the preferred embodiment allows for a substantial speed-up in the 
calculation of aggregate-layer transfer matrices [IV*//]. For instance, each parallel 
processor can be assigned to generate aggregate-layer transfer matrices UV^] over a 
separate range of wavelengths X since calculations at one wavelength X are independent 
of calculations at other wavelengths X, The I/O equipment 1505 can also include a 
keyboard 1502 and mouse 1504 for the input of information, and a display 1501 and 
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printer 1503 for the output of information. Many variations of this simple computer 
system 1500 can include systems with multiple I/O devices, multiple computers 
connected by Internet linkages, multiple computers connected by local area networks, 
devices that incorporate CPU 1515, memory 1520, and I/O equipment 1505 into a single 
unit, etc. 

[0082] As will be described below, in the present exemplary embodiment, 

portions of the above-described analysis are pre-computed and cached, thereby reducing 
the computation time required to calculate the diffracted reflectivity produced by a 
periodic grating. Briefly, the pre-computation and caching can include: 

• pre-computation and caching (i.e. , storage in a look-up 
table) of the permittivity s(x), the harmonic components Si of the permittivity s(x) 
and the permittivity harmonics matrix [E], and/or the permittivity s(x), the 
inverse-permittivity harmonics % t and the inverse-permittivity harmonics matrix 
[P] 9 for intra-layer parameters spanning an intra-layer parameters region 

• pre-computation and caching of the wave-vector matrix [A] 
for intra-layer parameters spanning an intra-layer parameters region {ju} and 
incident-radiation parameters spanning an incident-radiation parameters region 

{k}; 

• pre-computation and caching of eigenvectors w i>m and 
eigenvalues x m of the wave-vector matrix [A] to form an eigenvector matrix [W\ 9 a 
root-eigenvalue matrix [Q] 9 and a compound matrix [V] for intra-layer parameters 
spanning an intra-layer parameters region {\x} and incident-radiation parameters 
spanning an incident-radiation parameters region {k}; and 

• pre-computation and caching of aggregate-layer transfer 
matrices [r^y] for aggregate-layer parameters spanning an aggregate-layer 
parameters region {<)>}. 
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[0083] Briefly, in the present exemplary embodiment, the use of the pre- 

compiled and cached aggregate-layer transfer matrices [IV^] to calculate the diffraction 
spectrum from a periodic grating can include: 

• construction of matrix equation (2. 1 .9) or (2. 1 . 1 0) by 
retrieval of cached aggregate-layer transfer matrices [r^j] corresponding to the 
aggregate-layer trapezoids 1030 of the grating 145 under consideration; and 

• solution of the matrix equation (2. 1 .9) or (2. 1 . 1 0) to 
determine the diffracted reflectivity R( for each harmonic order /. 

[0084] The intra-layer precalculation and caching portion is illustrated by 

consideration of the exemplary ridge profiles 1601 and 1651 shown in cross-section in 
Figs. 16A and 16B, respectively. The profile 1601 of Fig. 16A is approximated by four 
slabs 1611, 1612, 1613 and 1614 of rectangular cross-section. Similarly, the profile 1651 
of Fig. 16B is approximated by three slabs 1661, 1662, and 1663 of rectangular cross- 
section. The two exemplary ridge profiles 1601 and 1651 are each part of an exemplary 
periodic grating (other ridges not shown) which have the same grating period D, angle 
9 of incidence of the radiation 910, incident radiation wavelength X 9 and the index of 
refraction n 0 of the atmospheric material between the ridges 1601 and 1651 is the same. 
It should be noted that slabs 1614 and 1663 do not have the same thicknesses t, nor do 
slabs 1613 and 1661, or slabs 1611 and 1662 have the same thicknesses t However, 
slabs 1613 and 1661 have the same ridge slab width d 9 x-offset p, and index of refraction 
n r . Similarly, slabs 1611 and 1662 have the same ridge slab width d, x-offset p, and 
index of refraction n r , and slabs 1614 and 1663 have the same ridge slab width d 9 x-offset 
P, and index of refraction n r . Note that the ridges can be mounted directly on a substrate, 
or on films deposited on a substrate, since a film can be considered to be a ridge having a 
width d equal to the pitch D. 

[0085] Note that thickness t is not a parameter upon which the wave-vector 

matrix [A] is dependent. Although thickness t is actually an intra-layer parameter, for 
reasons that will become clear below, thickness t will be considered an aggregate-layer 
parameter. Therefore, the wave-vector matrices [A] are the same for slabs 1611 and 
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1662, 1613 and 1661, and 1614 and 1663, and the eigenvector matrices [FT], the root- 
eigenvalue matrices [Q] 9 and the compound eigensystem matrices [V\ for slabs 1661, 
1662 and 1663 are the same as the eigenvector matrices [W], the root-eigenvalue matrices 
[Q], and the compound eigensystem matrices [V\ for slabs 1613, 161 1 and 1614, 
respectively. Therefore, caching and retrieval of the eigensystem matrices [W], [Q], and 
[V\ for slabs 1613, 1611 and 1614 prevents the need for recalculation of eigensystem 
matrices [W], [Q], and [V] for slabs 1661, 1662 and 1663, and can reduce the 
computation time. Additionally, the pre-calculation and caching of eigensystem matrices 
[£?]> an d W\ f° r useful ranges and samplings of intra-layer parameters and incident- 
radiation parameters can reduce the computation time. 

[0086] As can be seen from equations (1.1.2), (1.1.3), (2.1.2) and (2.1.3), the 

permittivity harmonics s^- and the inverse permittivity harmonics tzij are only dependent 
on the following intra-layer parameters: the index of refraction of the ridges n Ti the index 
of refraction of the atmospheric material n 0 , the pitch Z>, the ridge slab width d, and the jc- 
offset p. As shown in the flowchart of Fig. 13, method (1300) begins with the 
determination (1305) of the ranges n rMn to n Kmax , n 0Mn to n 0>max , D min to D max , d min to d max , 
and $ min to Pmax, and sampling increments hn n 5n 0 , 5D, bd, and 8p for the intra-layer 
parameters, i.e., the index of refraction of the ridges n r the index of refraction of the 
atmospheric material n 0 , the pitch D, the ridge slab width d y and the x-offset p. The 
determination (1305) of the maximum harmonic order o must also be made. This 
information is forwarded from the I/O device 1505 to the CPU 1515. Typically, when 
applied to periodic gratings produced by semiconductor fabrication techniques, the 
ranges n r)min to n r , m ax, n 0f mm to n 0 ,max, D min to D max , d min to d max , and $ min to $ max are 
determined based on knowledge and expectations regarding the fabrication materials, the 
fabrication process parameters, and other measurements taken of the periodic grating 145 
or related structures. Similarly, when matching a measured diffraction spectrum to a 
library of calculated diffraction spectra to determine the dimensions of the periodic 
grating that created the measured diffraction spectrum, the sampling increments 8w r , 5n 0 , 
5D, hd, and 5p, and maximum harmonic order o of the library, are chosen based on the 
resolution to which the intra-layer parameters n r> no, A d and p are to be determined. The 
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intra-layer parameter ranges n r , min to n rjnax , n 0Mn to n 0 , ma x, D min to D max , d min to d max , and 
Pmrn to $ mox , and increments bn r , 8«o, 5D, 6d, and 5p define a five-dimensional intra-layer 
caching array More specifically, the caching array {jj,} consists of intra-layer points 
with the n r coordinates being {n r , m i n , n r>min +bn r , n nm i n +2&n ri n rtmax -25n rj n r>max -§n r> 
n^max}, the n 0 coordinates being {n 0ymini n 0 , m m+§no, n 0 , m w+2dno, n 0 , m ax-28n 0 , no,max-Sn 0 , 
n 0fmax }, the D coordinates being {D mim A™+8£>, D min +2§D, . . D max -2hD t D max -§D D max }, 
the d coordinates being {d mim d m i n +bd t d min +28d f d max -2hd j d max -hd t d max ), and the p 
coordinates being {$ min , P^+dp, Pm/«+25p, p mflX -25p, p W o*-5p, p wax }. In other words, 
the intra-layer caching array {ji} is defined as a union of five-dimensional coordinates as 
follows: 

{/A = U (w r min + i*i r , w 0>min + jSn 0 , + kSD, d^+tnSdtjS^+ndfi), 

i,j,k,l,m 

(3.1.1) 

where i,j, k, m and « are integers with value ranges of 

0<i<(n r ^ x -n rmiB )/5n r , (3.1.2a) 

0<j<(n 0 ^-n OM Jidn 0 , (3.1.2b) 

OZkZiD^-D^/dD , (3.1.2c) 

0<m<{d wx -d^)IM, (3.1.2d) 

and 

OZnZifi^-fi^yW. (3.1.2.e) 

It should be noted that the layer subscript, /, is not used in describing the intra-layer 
parameters n ry no, D 9 d, and p used in the layer-property caching array {\x} because each 
particular point juj in the layer-property caching array may correspond to none, one, 
more than one, or even all of the layers 1025 of a particular periodic grating 145. 



35 



[0087] As shown in Fig. 13, for each point in the layer-property caching array 

the "formulation-required" permittivity harmonics z i die calculated (1310) by the 

CPU 1515 and cached (1315) in memory 1520, and the "formulation-required" 

permittivity harmonics matrices p^are compiled from the cached formulation-required 

permittivity harmonics e . and cached (1315') in memory 1520. For RCWA analyses of 

TE-polarized incident radiation 910, or RCWA analyses of TM-polarized incident 
radiation 910 according to the formulation of equations (1.5.3) and (1.5.1 1), the 

formulation-required permittivity harmonics s i are the permittivity harmonics e, 
calculated (1310) according to equations (1.1.2) and (1.1.3), and the formulation-required 
permittivity harmonics matrix \e] is the permittivity harmonics matrix [E] formed as per 
equation (1.1.4). 

[0088] Similarly, for RCWA analyses of TM-polarized incident radiation 910 

according to the formulation of equations (1.5.1) and (1.5.9) or equations (1.5.2) and 

( 1 . 5 . 1 0), the formulation-required permittivity harmonics s f are the permittivity 

harmonics 8/ calculated (1310) according to equations (1.1.2) and (1.1.3) and the inverse- 
permittivity harmonics 71/ calculated (1310) according to equations (1.5.6) and (1.5.7), 
and the formulation-required permittivity harmonics matrices \e] are the permittivity 
harmonics matrix [E] formed from the permittivity harmonics e* as per equation (1.1.4) 
and the inverse-permittivity harmonics matrix [P] formed from the inverse-permittivity 
harmonics 71/ as per equations (1.5.6) and (1.5.7). 

[0089] As per equations (1.3.7), (1.5.1), (1.5.2) and (1.5.3), the wave-vector 

matrix [A] is dependent on the required permittivity harmonics matrices \e\ and the 
matrix [K x ]. The matrix [K x ], in addition to being dependent on "global" layer-property 
parameters (i.e., the atmospheric index of refraction n 0 and pitch D), is dependent on 
incident-radiation parameters, i.e., the angle of incidence 0 and the wavelength X of the 
incident radiation 910. As shown in the flowchart of Fig. 13, ranges 0 m i n to 0 max and X m \ n 
to A,max, and increments 50 and 5 A, are determined (1317) for the incidence angle 0 and 
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wavelength X, and forwarded from the I/O device 1505 to the CPU 1515. The incident- 
radiation caching array {k} is defined as a union of two-dimensional coordinates as 
follows: 

{*} = U +pS0, 1^ + qhX) (3.1.3) 

n,o 

where p and q are integers with value ranges of 

OZpZ(0 m -O^)/8B , (3.1.4a) 

0<q<{X^-X^)lbX. (3.1.4b) 

[0090] The combined caching array {p., K}is formed (1307) by the union of 

coordinates as follows: 

{//, k) = U (n rMn + iSn r , « 0>mill + y7Si 0 , Z) min + MD, rf^ + m&t, y? min + 0 min + /><?0, + qdX) 

i,j,k,l,m 

(3.1.5) 

where i, j, k, m, n,pzndq satisfy equations (3.1.2a), (3.1.2b), (3.1.2c), (3,1. 2d), (3.1.4a) 
and (3.1 .4b). Typically, the ranges 0 min to 9 max and A, min to X max are determined (1317) 
based on knowledge and expectations regarding the apparatus (not shown) for generation 
of the incident radiation 910 and the apparatus (not shown) for measurement of the 
diffracted radiation 914. Similarly, the sampling increments 50 and 8 A, are determined 
(1317) based on the resolution to which the layer-property parameters n r> n 0 , D, d, and P 
are to be determined, and/or the resolution to which the incident-radiation parameters 9 
and X can be determined. For each point in the combined caching array k}, the wave- 
vector matrix [A] is calculated (1320) by the CPU 1515 according to equation (1.3.7), 
(1.5.1), (1.5.2) or (1.5.3) and cached (1225). 

[0091] Since the wave-matrix matrix [A{\ is only dependent on incident-radiation 

parameters (angle of incidence 0 of the incident radiation 910, wavelength X of the 
incident radiation 910) within layer /, and intra-layer parameters (index of refraction of 
the ridges n r> index of refraction of the atmospheric material n 0 , pitch D, ridge slab width 
d, x-offset (5), and it follows that the eigenvector matrix [W{\ and the root-eigenvalue 
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matrix [Qi\ are also only dependent on the incident-radiation parameters 9 and X, and the 
layer-property parameters n r , no, A d, and (5 within layer /. According to a preferred 
embodiment, the eigenvector matrix [W{\ and its root-eigenvalue matrix [Qi] are 
calculated (1347) by the CPU 1515 and cached (1348) in memory 1520 for each point in 
the combined caching array {\x y k}. The calculation (1347) of the eigenvector matrices 
[W{] and the root-eigenvalue matrices [Qi\ can be performed by the CPU 1515 using a 
standard eigensystem solution method, such as singular value decomposition (see 
Chapter 2 of Numerical Recipes, W.H. Press, B.P. Glannery, S.A. Teukolsky and W.T. 
Vetterling, Cambridge University Press, 1986). The matrix [V], where [V]= [W] [Q], is 
then calculated (1357) by the CPU 1515 and cached (1358) in memory 1520. 

[0092] Then, a library of useful aggregates of layers is generated. It should be 

noted that although the use of aggregate-layer trapezoids has been described, useful 
aggregates of layers may have other shapes as well. For instance, useful aggregates of 
layers may have more than four corners, may have convex as well as concave corners, 
may have curved sides, etc. Generally, aggregate-layer trapezoids will have horizontal 
top and bottom edges, so that boundaries between aggregate-layer trapezoids will 
coincide with boundaries between harmonic expansion layers 1025. 

[0093] As shown in Fig. 17, aggregate-layer trapezoids are described by four 

aggregate-layer parameters: the width b of the base 1705 of the trapezoid 1030, the width 
w of the top 1710 of the trapezoid 1030, the height h of the trapezoid 1030, and the 
distance offset c between the center 1707 of the base 1705 of the trapezoid 1030 and 
center 1712 of the top 1710 of the trapezoid 1030. Alternatively, the offset between the 
center 1707 of the base 1705 of the trapezoid 1030 and center 1712 of the top 1710 of the 
trapezoid 1030 may be given by a parameter which is an angle or a dimensionless ratio. 

[0094] The left ridges of Figs. 10B and 10C show groupings of the rectangular 

slabs 1026 into aggregate-layer trapezoids 1030 with side edges which approximately 
follow the side contours of the rectangular slabs 1026. In Fig. 10B, rectangular slabs 
1026.1 and 1026.2 are grouped into the upper aggregate-layer trapezoid 1030.1, 
rectangular slabs 1026.3 and 1026.4 are grouped into the center aggregate-layer trapezoid 
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1030.2, which happens to be rectangular, and rectangular slabs 1026.5, 1026.6 and 
1026.7 are grouped into the lower aggregate-layer trapezoid 1030.3. It should be noted 
that the upper and lower trapezoids 1030.1 and 1030.3 have side edges that are not 
vertical, while the rectangular slabs 1026 in each group represented by a trapezoid 1030.1 
and 1030.3 have vertical side edges. The number of aggregate-layer trapezoids into 
which the rectangular slabs 1026 are grouped, and the particular assignment of 
rectangular slabs 1026 to aggregate-layer trapezoid 1030.3, will depend on the algorithm 
used to determine the groupings. For instance, in Fig. 10C, rectangular slab 1026.1 is 
represented by the upper, rectangular aggregate-layer trapezoid 1030.1, rectangular slabs 
1026.2 and 1026.3 are grouped into the center aggregate-layer trapezoid 1030.2, and 
rectangular slabs 1026.4, 1026.5, 1026.6 and 1026.7 are grouped into the lower 
aggregate-layer trapezoid 1030.3. In a very crude approximation, rectangular slabs 
1026.1 through 1026.7 could be represented by a single aggregate-layer trapezoid. 
Generally, the accuracy of a representation with aggregate-layer trapezoids will increase 
with the number of aggregate-layer trapezoids that are used, until the number of 
aggregate-layer trapezoids approaches the number (L-l) of rectangular slabs 1026 in a 
ridge 902. Typically, while between 16 and 25 rectangular slabs 1026 are required to 
accurately represent a profile, only 2 to 4 aggregate-layer trapezoids 1030 are required to 
accurately represent a profile. Therefore, even ignoring the speed-up due to caching, the 
calculation of the matrix product according to equations (2.1.9) or (2.1.10) for three 
aggregate-layer transfer matrices can be seven times faster than the calculation of 

the matrix product according to equations (2.1.5) or (2.1.6) for 21 non-aggregate layer 
transfer matrices |T/]. 

[0095] The aggregate-layer precalculation and caching portion is illustrated by 

consideration of the exemplary ridge profiles 1678 and 1679 shown in cross-section in 
Figs. 16C and 16D, respectively. The profile 1678 of Fig. 16A is approximated by four 
aggregate-layer trapezoids 1671, 1672, 1673, and 1674. Similarly, the profile 1679 of 
Fig. 16D is approximated by four aggregate-layer trapezoids 1671, 1672, 1673, and 1676. 
The two exemplary ridge profiles 1678 and 1679 are each part of exemplary periodic 
gratings (other ridges not shown) which have the same grating period D, angle 9 of 
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incidence of the radiation, incident radiation wavelength X 9 and index of refraction n 0 of 
the atmospheric material. The aggregate-layer trapezoid 1671 has the same base width b, 
top width w, height h and top-base offset c in both Figs. 16C and 16D. Furthermore, 
although not discernible from Figs. 16C and 16D, the index of refraction n r of aggregate- 
layer trapezoid 1671 is the same in both Figs. 16C and 16D. Similarly, aggregate-layer 
trapezoids 1672 and 1673 have the same base width b, top width w, height top-base 
offset c, and index of refraction n r in both Figs. 16C and 16D. However, aggregate-layer 
trapezoids 1674 and 1676 of Figs. 16C and 16D have no equivalents in Figs. 16D and 
16C, respectively. 

[0096] If aggregate-layer trapezoid 1671 in Fig. 16C and aggregate-layer 

trapezoid 1671 in Fig. 16D are composed of sets of rectangular slabs of the same 
dimensions (i.e., the same thicknesses t and widths d) which are stacked in the same 
order, then, as per equation (2.1.7), the aggregate-layer transfer matrix [rV»/y] is the same 
for aggregate-layer trapezoid 1671 in Fig. 16C and aggregate-layer trapezoid 1671 in Fig. 
16D. Similarly, if aggregate-layer trapezoids 1672 and 1673 in Fig. 16C and aggregate- 
layer trapezoids 1672 and 1673 in Fig. 16D are composed of sets of rectangular slabs of 
the same dimensions which are stacked in the same order, then the aggregate-layer 
transfer matrix [Tu^y] is the same for the trapezoids 1672 and 1673 in Fig. 16C and the 
trapezoids 1672 and 1673 in Fig. 16D, respectively. Therefore, caching and retrieval of 
the aggregate-layer transfer matrices [Tn^y] for aggregate-layer trapezoids 1671, 1672 
and 1673 prevents the need for recalculation of the aggregate-layer transfer matrices 
[Tii^ij] and reduces the computation time. Additionally, the pre-calculation and caching 
of aggregate-layer transfer matrices [Tu^y] for useful ranges and samplings of incident- 
radiation parameters, intra-layer parameters and aggregate-layer parameters can 
substantially reduce the computation time when multiple spectra utilizing those 
aggregate-layer transfer matrices [T^y] are to be generated. 

[0097] The mapping from rectangular slabs to aggregate-layer trapezoids 1030 is 

accomplished by a least-squares fit of the sides of the rectangular slabs to straight lines. 
Initially, the lengths of the fit lines are a substantial fraction of the total height of the 
profile. However, if the goodness-of-fit exceeds a critical value, the length of a fit line is 
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reduced until the goodness-of-fit is below the critical value. It should be noted that the 
mapping from sets of stacked rectangular slabs to aggregate-layer transfer matrices 
[Tii^ij] is generally many-to-one. For instance, for the rectangular trapezoid 1672, the 
rectangular slabs that it represents might all have sides which are well-aligned with the 
sides of the rectangular trapezoid 1672. Or, the left and right sides of the rectangular 
slabs might deviate somewhat from the left and right sides of the rectangular trapezoid 
1672, but have an average left-side position coincident with the left side of the trapezoid 
1672 and an average right-side position coincident with the right side of the trapezoid 
1672. Therefore, due to the many-to-one mapping, a retrieved aggregate-layer transfer 
matrix [IV^] generally only provides an approximation to the aggregate-layer transfer 
matrix [TV>/y] that would actually be generated by the set of rectangular slabs . 

[0098] As shown in the flowchart of Fig. 1 3, the ranges w min to w maX9 b min to b max , 

hmin to h max , and c min to c max , and sampling increments 5w, 56, 8A, and 5c for the 
aggregate-layer parameters {i.e., the base width 6, the top width w, the height h, and the 
base-top offset c) are determined (1360), and forwarded from the I/O device 1505 to the 
CPU 1515. Typically, when applied to periodic gratings produced by semiconductor 
fabrication techniques, the ranges w min to w max , b m in to b maXy h min to h max , and c min to c max 
are determined based on knowledge and expectations regarding the fabrication materials, 
the fabrication process parameters, and other measurements taken of the periodic grating 
145 or related structures. Similarly, when matching a measured diffraction spectrum to a 
library of calculated diffraction spectra to determine the dimensions of the periodic 
grating that created the measured diffraction spectrum, the sampling increments 8w 56, 
hh, and 8c are chosen based on the resolution to which the aggregate-layer parameters w, 
b y h and c are to be determined. The aggregate-layer parameter ranges w min to w maX9 b min 
to b max , h min to h maX9 and c min to c max , and increments 5w, 56, 5/z, and 8c define a four- 
dimensional aggregate-layer caching array {<))}. More specifically, the aggregate-layer 
caching array {(()} consists of aggregate-layer points with the w coordinates being {w min , 
w OT ^+5w, w min +2$w, . . . 3 Wmax-§w, w max }, the b coordinates being {b miny b min +&b, b min +2bb, 
. . . , bmax~§b f b max ) , the h coordinates being {h min , h min +dh, h min +2bh, . . . , h max -2bh, h max - 
5/2, h max }, and thee coordinates being {c min , c min +8c, c min +28c c max -28c, c max -Sc, c max } . 
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In other words, the aggregate-layer caching array {<[>} is defined as a union of four- 
dimensional coordinates as follows: 




(3.1.6) 



where i, j, k and m are integers with value ranges of 



0<i<(w T 



max 



-w in )/5w , 



(3.1.7a) 



o<y<(^ 



(3.1.7b) 



0<k<(h, 



max 



(3.1.7c) 



and 



0<m<(c, 



max 



(3.1.7 A) 



[0099] 



It should be noted that if any of the intra-layer parameters n r , no, D 9 d, and 



p, any of the incident-radiation parameters 9 and X, or any of the aggregate-layer 
parameters b, w, h and c are known to sufficient accuracy, then a single value, rather than 
a range of values, of the parameter may be used, and the dimensionality of the 
corresponding caching array {p,}, {k} or {§} is effectively reduced. It should also be 
understood that the incident-radiation parameters region, the intra-layer parameters 
region and/or the aggregate-layer parameters region need not be hyper-rectangles, and the 
incident-radiation parameters region, the intra-layer parameters region and/or the 
aggregate-layer parameters region need not be sampled on a uniform grid. For instance, a 
sampling may be performed using a stochastic sampling method. Furthermore, the 
sampling density need not be uniform. For instance, the sampling density may decrease 
near the boundaries of a parameters region if circumstances near the boundaries are less 
likely to occur. 

[001 00] Once the range and sampling for the aggregate-layer caching array {$} is 
determined (1360), layer transfer matrices [F{\ are calculated (1370) by the CPU 1515 for 
each aggregate-layer trapezoid 1030, and cached (1375) in memory 1520 in the 
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aggregate-layer caching array {<)>} . As shown in Fig. 13, the calculation of the layer 
transfer matrices [T{\ utilizes the cached eigenvectors [W{\, root-eigenvalues [£)/], and 
compound eigensystem matrices [Vfl. Furthermore, due to the dependence of the layer 
transfer matrix [T/] on the layer translation matrix [X/], as per equation (2.1 .4), the 
thickness / of the rectangular slab 1026 is provide for the calculation. As mentioned 
above, although thickness t is actually an mtra-layer parameter, thickness t is considered 
an aggregate-layer parameter for the purposes of the present specification since it is not 
utilized until the calculation (1370) of the layer transfer matrices |T/]. Utilizing the 
cached layer transfer matrices |T/], aggregate-layer transfer matrices [Tu^iJ] are 
calculated (1370) by the CPU 1515 for each point in the aggregate-layer caching array 
{(|)}, and cached (1385) in memory 1520. 

[00101] The use of pre-computed and cached aggregate-layer transfer matrices 
[Tif^ij] is shown in Fig. 14. Use of the cached aggregate-layer transfer matrices [TV^] 
begins by aggregate-layer discretizing a ridge profile and determining the parameters 
describing the aggregate-layer discretized ridge profile. In particular, the incident- 
radiation parameters (i.e., the angle of incidence 0 and the wavelength X of the incident 
radiation), the intra-layer parameters (i.e., index of refraction of the ridges « r> the index of 
refraction of the atmospheric material n 0 , the pitch Z>, the ridge slab width d 9 and the x- 
offset P) for each rectangular slab 1026 in each aggregate-layer trapezoid 1030, and the 
aggregate-layer parameters (i.e., the base width b, the top width w 9 the height h, and the 
base-top distance offset c) for each aggregate-layer trapezoid 1030 are determined (1405) 
and forwarded via the I/O device 1505 to the CPU 1515. 

[001 02] Once the incident-radiation, intra-layer, and aggregate-layer parameters 
for an aggregate-layer discretized profile are determined (1405), the cached aggregate- 
layer transfer matrices |JV»/y] for the aggregate-layer trapezoids of the discretized ridge 
profile are retrieved (1410) from memory 1520 for use by the CPU 1515 in constructing 
(1415) the system matrix equation (2.1.9) or (2.1.10). The CPU 1515 then solves (1420) 
the boundary-matched system matrix equation (2.1.9) or (2.1.10) for the reflectivity i? z of 
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each harmonic order from -o to +o and each wavelength X of interest, and forwards the 
results to an output device 1505 such as the display 1501, printer 1503, or the like. 



[001 03] Fig. 1 8 presents a flowchart for the use of a library of profile-spectrum 
pairs, and use of the library to determine the profile corresponding to a measured 
spectrum. As shown in Fig. 18, the numerical aperture of the focusing optics in the 
optical path of the reflectometry or ellipsometry configuration is characterized (1810), 
and the range and sampling of incident-radiation parameters {i.e., incidence angle 0 and 
wavelength X) required for the calculation of the simulated diffraction spectra are 
determined (1812). Ranges and samplings of intra-layer parameters {i.e., the index of 
refraction of the ridges w r> the index of refraction of the atmospheric material n 0 , the pitch 
A the ridge slab width d, and the x-offset |3) corresponding to expected ranges of profile 
shapes are also determined (1815). The intra-layer dependent portions of the rigorous 
coupled-wave calculation {i.e., the wave-vector matrices [A], eigenvector matrices [W\ 9 
root-eigenvalue matrices [Q] 9 and compound matrices [V]) are then pre-computed (1820) 
by the CPU 1515 for the determined ranges and samplings of the intra-layer parameters 
and incident-radiation parameters, and cached (1825) in the cache memory 1520. Ranges 
and samplings of aggregate-layer parameters {i.e., base width b, top width w, height h, 
and base-top offset c) corresponding to expected ranges of profile shapes are also 
determined (1827). Then, using the pre-computed and cached intra-layer portions of the 
rigorous coupled-wave calculation, aggregate-layer transfer matrices [T/^/y] 
corresponding to the determined range and sampling of aggregate-layer parameters are 
calculated (1828) by the CPU 1515 and cached (1829) in the cache memory 1520. 
Diffraction spectra are then calculated (1830) for combinations of aggregate-layer 
trapezoids 1030 to produce a sub-library of spectra. Each combination of aggregate-layer 
trapezoids 1030 is indexed (1835) in the cache memory 1520 with its associated spectrum 
to provide a library of profile-spectrum pairs. Diffraction signals are acquired (1 805) by 
an optical setup, using either a reflectometry or an ellipsometry configuration, as depicted 
in Figs. 1 and 2. The diffraction signal may be acquired prior to, during, or subsequent to 
the generation (1810), (1820), (1825), (1827), (1828), (1829), (1830) and (1835) of the 
profile-spectra library. Finally, a measured diffraction spectrum obtained from a physical 
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profile is compared (1840) to the spectra in the library in order to select the spectrum in 
the library that best matches the measured diffraction spectrum. The combination of 
aggregate-layer trapezoids 1030 corresponding to the best-matched calculated spectrum 
is the estimate of the physical profile. The matching algorithms that can be used for this 
purpose range from simple least-squares approach (linear regression) to a neural network 
approach that compares features of the measured signal with the spectra from the library 
using a non-linear, principal-component regression scheme. Explanations of such 
methods may be found in numerous text books on the topic, such as Chapter 14 of 
"Mathematical Statistics and Data Analysis" by John Rice, Duxbury Press and Chapter 4 
of "Neural Networks for Pattern Recognition" by Christopher Bishop, Oxford University 
Press. 

[001 04] It should be understood that the present invention is also applicable to off- 
axis or conical incident radiation 910 (i.e., the case where (j)^0 and the plane of incidence 
912 is not aligned with the direction of periodicity, x , of the grating). The above 
exposition is straightforwardly adapted to the off-axis case since, as can be seen in 
"Rigorous Coupled- Wave Analysis of Planar-Grating Diffraction," M. G. Moharam and 
T. K. Gaylord, J. Opt. Soc. Am., vol. 71, 81 1-818, July 1981, the differential equations 
for the electromagnetic fields in each layer have homogeneous solutions with coefficients 
and factors that are only dependent on intra-layer parameters and incident-radiation 
parameters. As with the case of on-axis incidence, aggregate-layer calculations are pre- 
calculated and cached. In computing the diffracted reflectivity from a periodic grating, 
cached aggregate-layer transfer matrices [rV»/ y ] corresponding to aggregate-layer 
trapezoids of the profile are retrieved for use in constructing a system matrix equation in 
a manner analogous to that described above. 

[001 05] Although the use of rigorous coupled-wave method of calculating the 
diffraction of radiation has been described, any method of diffraction calculation where 
the system equation involves the multiplication of a series of matrices, where each of 
these matrices is dependent on local parameters (such as intra-layer parameters) and/or 
global parameters (such as incident-radiation parameters) can be used. For example, the 
diffraction calculation may be an approximate method, and/or an integral formulations, or 
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any other formulations, such as those mentioned in standard texts such as Solid State 
Physics, N.W. Ashcroft and ND. Mermin, Saunders College, Philadelphia, 1976, pages 
133-134, or Optical Properties of Thin Solid Films, O.S. Heavens, Dover Publications, 
Inc., New York, 1991, or Ellipsometry and Polarized Light, R.M.A. Azzam and N.M. 
Bashara, North-Holland Personal Library, Amsterdam, 1987. Furthermore, the present 
invention maybe applied to diffraction calculations based on decompositions or analyses 
other than Fourier analysis, such as a decomposition into Bessel functions, Legendre 
polynomials, wavelets, etc. 

[001 06] The foregoing descriptions of specific embodiments of the present 
invention have been presented for purposes of illustration and description. They are not 
intended to be exhaustive or to limit the invention to the precise forms disclosed, and it 
should be understood that many modifications and variations are possible in light of the 
above teaching. The embodiments were chosen and described in order to best explain the 
principles of the invention and its practical application, to thereby enable others skilled in 
the art to best utilize the invention and various embodiments with various modifications 
as are suited to the particular use contemplated. Many other variations are also to be 
considered within the scope of the present invention. 

[001 07] For instance: aggregate-layer transfer matrices may be generated and 
cached without caching of intra-layer calculations; the calculation of the present 
specification is applicable to circumstances involving conductive materials; once the 
aggregate-layer transfer matrices are calculated and cached, intermediate results (such as 
the permittivity, inverse permittivity, permittivity harmonics, inverse-permittivity 
harmonics, permittivity harmonics matrix, the inverse-permittivity harmonics matrix, the 
wave-vector matrix, eigenvalues of the wave- vector matrix, and/or eigenvalues of the 
wave-vector matrix) need not be stored; the compound matrix [FJ, which is equal to the 
product of the eigenvector matrix and the root-eigenvalue matrix, may be calculated 
when it is needed, rather than cached; the eigenvectors and eigenvalues of the matrix [A] 
may be calculated using another technique; a range of an incident-radiation parameter, 
intra-layer parameter and/or aggregate-layer parameter may consist of only a single 
value; the array of regularly- spaced incident-radiation parameter values, intra-layer 
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parameter values and/or aggregate-layer parameter values for which the matrices, 
eigenvalues and eigenvectors are cached may be replaced with an array of irregularly- 
spaced incident-radiation parameter values, intra-layer parameter values and/or 
aggregate-layer parameter values, or a random selection of incident-radiation parameter 
values, intra-layer parameter values and/or aggregate-layer parameter values; the system 
equation may be solved for the diffracted reflectivity and/or the diffracted transmittance 
using any of a variety of matrix solution techniques; the "ridges" and "troughs" of the 
periodic grating may be ill-defined; the method of the present invention may be applied 
to gratings having two-dimensional periodicity; the method of the present invention may 
be applied to any polarization which is a superposition of TE and TM polarizations; the 
ridged structure of the periodic grating may be mounted on one or more layers of films 
deposited on the substrate; the method of the present invention may be used for 
diffractive analysis of lithographic masks or reticules; the method of the present 
invention may be applied to sound incident on a periodic grating; the method of the 
present invention may be applied to medical imaging techniques using incident sound or 
electromagnetic waves; the method of the present invention may be applied to assist in 
real-time tracking of fabrication processes; the gratings may be made by ruling, blazing 
or etching; the grating may be periodic on a curved surface, such as a spherical surface 
or a cylindrical surface, in which case expansions other than Fourier expansions would be 
used; the method of the present invention may be utilized in the field of optical analog 
computing, volume holographic gratings, holographic neural networks, holographic data 
storage, holographic lithography, Zernike's phase contrast method of observation of 
phase changes, the Schlieren method of observation of phase changes, the central dark- 
background method of observation, spatial light modulators, acousto-optic cells, etc. In 
summary, it is intended that the scope of the present invention be defined by the claims 
appended hereto and their equivalents. 
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